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Abstract 

Two different local stability conditions for an asymmetric superfluid has been discussed in the 
literature. We here consider the relations between them. 
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In our recent paper [1], we considered a two component Fermi gas with unequal popula- 
tions under a short-ranged attractive interaction. We constructed the phase diagram by first 
: inding uniform state solutions to the number and gap equations. We then impose the local 
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stability condition that the susceptibility matrix has to have only positive eigenvalues. 
In particular then, in terms of the differences in populations and chemical potentials rid and 
h, local stability /instability is decided via 

dnA > (1) 



dh < 

Here /i is the average chemical potential, [.3] There we stated, without providing details, 
that a negative susceptibility implies instability towards phase separation. We thus then 
exclude all states with negative susceptibilities and identified the corresponding point in the 
phase diagram as not having a stable uniform state. 

The phase diagram and local stability condition have been discussed in a different manner 
flflflflflafi These authors considered the fr ee e„e rgy of the 8y8 te,n fl(A; „, H) as 
a function of the order parameter A for fixed \i and h [11]. The locally stable states are 
selected according to the condition that Q must be a relative minimum of Q as a function 
of A. Thus, at the extremum points = 0, they apply the condition 

(2) 

h 

to distinguish locally stable and unstable states. 

Since our condition eq (JTJ) was stated without providing detailed justifications, there are 
some confusions as to how the conditions ([IJ and (J2J) are related. Sheehy and Radzihovsky 
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suggest recently that states with positive susceptibilities can still correspond to a 
relative maximum of the free energy 0(A; h). Other groups 3, 13 , 14 ] have set out to 
show the relation between the conditions JI} and (J2J). These proofs, while correct, involve 
algebra which motivation may be a priori somewhat mysterious. The purpose of this note 
is to provide perhaps a physically more transparent picture on the relations between the 
conditions (JTJ and (J2J. 

We start by considering the locus of the relative extremum of f2(A,//, h), that is, the 
solution to the gap equation ff^H = 0- This gives A as a function of \i and h. For 
defmiteness [3], we consider a "cut" of this function at some given \i. First, let us note 



that the normal state A = is always a solution to this gap equation, and hence always 
an extremum of Q (when considered a function of A at fixed \x and h). For some region 
of parameters, there may be in addition non-trivial (A 7^ 0) solution to the gap equation. 
An example deep in the BEC (strong attractive interaction) limit is as sketched in the inset 
of Fig n The non-trivial solution to the gap equation gives a A that is monotonically 
decreasing with h. (This non-trivial solution for h <h\ represents a gapless superfluid jl|). 
This function however is not necessarily single valued. An example is as shown in the main 
part of Fig [D This happens in particular in the weak interaction limit, as was first shown by 
Sarma jj| (and continues to hold until the attractive interaction is sufficiently strong). 

There, the non-trivial solution for A, instead of decreasing monotonically with h, "bends 
around" at the point B. For h < h^, there is one non-trivial solution (labelled by S-A), 
whereas there are two non-trivial solutions for 1%a < h < hs (labelled by A-B and B-C). 

First we note a simple but useful fact. The free energy Q must continue to increase 
as A — ► 00, as A = 00 is not a solution to the gap equation, and Q cannot decrease 
indefinitely. Since the solutions to the gap equation already give all the extrema, it follows 
that the solution to the gap equation with the largest A always corresponds to a relative 
minimum of Q. In fact, we have the following statements: 

(1) If A = is the only solution to the gap equation, then it is the free energy minimum. 

(2) If there is only one non-trivial solution A = A 1; then it is the free energy minimum, 
and the normal state is a relative maximum. (In this case, starting from A = 0, Q first 
decreases, then reaches its minimum at Ai, then increases for all A > Ai.) This applies to 
the inset for h < hi, and h <Ha for the main Fig^ 

(3) If there are two non-trivial solutions Ai < A 2 , then A x must be a relative maximum, 
whereas both A2 and the normal state are relative minima. (In this case, starting from 
A = 0, Cl first increases, then reaches a relative maximum at A l5 decreases till A 2 , then 
increases again for all A > A 2 . 

Therefore, for the non-trivial solutions in Fig^ the branch S-A is the absolute minimum 
for Q, A-B is a relative minimum, B-C must be a relative maximum. For the trivial solution 
A = 0, it corresponds to a relative maximum for h < he = Ha and a relative minimum for 
h > he- (see also the last paragraph of this note). 

Now we consider as a function of h. For the situation in the inset of Fig Q there 
is one branch of rid for the normal state, and one for the superfluid state. A sketch is as 
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shown in the inset of Fig El where these two branches are denoted by the dashed and full 
lines respectively. Note that the full line must have positive slope as it is stable. Therefore 
there is one possible rid for each h for h > h\, but two for each h < hi. Now consider 
decreasing the attractive interaction (or increasing ii) so that the A-h plot evolves from the 
form in the inset of Fig ^ to that in the main part of that figure. Consider now what must 
happen to the rid-h plot. There continues to be one branch (dashed line in main part of Fig 
I2J) which corresponds to the normal state. For the non-trivial solution to the gap equation 
(S-A-B-C in Fig^), there is an S-A-B part in Fig El with positive slope corresponding to 
the locally stable S-A-B part in Fig^ However, since there is one rid for each solution to 
the gap equation, there must be two possible n^'s corresponding to the non-trivial solutions 
in the range he = < h < tig = hp. Thus if we follow the path S-A-B-C in Fig[Q the 
corresponding path in the rid versus h plot must also "turn around"" at point B, where the 
A versus h curve turns around. At this point, (^§f) necessarily diverges and changes sign 
at B, As a result, the branch BC in FigEJhas {^§f) < 0, corresponding to ( ) < in 
FiglU 



At this point let us make connection with the arguments given in 
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In 



essence, these arguments are based on the following equation or its variation: 



dn d \ / dnA ( aA*)^ , , 

dh ) ,. ~{dh L. />fi N 
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Here fid is a function of A, fx, h defined by the relation fid = — (f^) • Note that, at the 

solution to the gap equation, n d (i-i, h) = n<i(A, /i, h) since rid = — (§?) = — ( §r ) since 

(m) = ®- Suppose we follow the line S-A-B-C in FigQJ At B, one changes from a branch 

where [J^J > to < 0. Since (f^lj appears in the denominator of eq Q, 
must changes from positive to negative through ±oo near B Q. As we have already seen 
above, this behavior is a natural consequence of the fact that the curve rid versus h must 
turn around at point B in Fig|21 

As seen from eq (J3J) however, since (^f) A > always (c.f. 0), there seems no general 
principle requiring that (7^) < whenever (f^I) < 0. Similarly, in principle with 
further change of interaction parameters, some part of the BC curve in Fig |2] may evolve 
further so that it can acquire a positive slope. However, for this to happen, it can easily be 
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seen that one must have some region of rid which corresponds to at least three possible /i's, 
and with at least one point where (^^) vanishes. This does not seem to happen for the 
current system of asymmetric Fermi gas in mean-field theory. Moreover, there is no reason 
why the change to positive slope would occur at the same point in parameter space as where 
A versus h turns around in Fig ^ Thus evaluating the susceptibility should still indicate 
when one goes from a local minimum to maximum of Q as a function of A. 

Finally, we notice here also that Fig |21 is similiar to the one in theory of liquid gas 
transition such as the Van der Waals equation of state for a non-ideal gas. h and — rid 
here plays the role of the pressure P and volume V respectively. (A sign change in rid is 
necessary since here dfl = —riddh, whereas for the gas, the differential of Gibbs free energy 
is dG = VdP). Indeed, Fig has precisely that form if we make the above identification 
(and ignore the branch NC, which is never stable). The branch BC there is locally unstable 
since the compressibility is negative. By again the analogy to the Van der Waals gas, the 
absolutely stable branch can be determined by Maxwell construction. These branches are 
S-A-E s for h < He = h,Es = and E n -D-N' if h > Iie, with the condition that the areas 
E s -B-M and M-C-E n are equal. At h = He, the system phase separates into a mixture with 
part of the volume in state E s and part in E n . 

This research was supported by the NSC of Taiwan under grant number NSC94-2112-M- 
001-002. 
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FIG. 1: A versus h at a given chemical potential. 
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FIG. 2: rirf versus h at a given chemical potential. 



